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Functional equations like exact renormalization group and Dyson-Schwinger equations have con- 
tributed to a better understanding of non-perturbative phenomena in quantum field theories in 
terms of the underlying Green functions. In Yang-Mills theory especially the Landau gauge has 
been used, as it is the most accessible gauge for these methods. The growing understanding ob- 
tained in this gauge allows to proceed to other gauges in order to obtain more information about 
the relation of different realizations of the confinement mechanism. In the maximally Abelian 
gauge first results are very encouraging as a variant of Abelian infrared dominance is found: The 
Abelian part of the gauge field propagator is enhanced at low momenta and thereby dominates 
the dynamics in the infrared. Its role is therefore similar to that of the ghost propagator in the 
Landau gauge, where one denotes the corresponding phenomenon as ghost dominance. Also the 
ambiguity of two different types of solutions (decoupling and scaling) exists in both gauges. Here 
we present how the two solutions are related in the maximally Abelian gauge. The intricacy of 
the system of functional equations in this gauge required the development of some new tools and 
methods as, for example, the automated derivation of the equations by the program DoFun. We 
also present results for linear covariant and ghost anti-ghost symmetric gauges. 
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1. Dual superconductor and Abelian infrared dominance 

Despite the success of quantum chromodynamics in describing many aspects of the strong 
interaction we still lack a complete understanding of the phenomenon of confinement in terms 
of the underlying fundamental fields, the quarks and gluons. Many scenarios have been brought 
forward, among these the dual superconductor picture [1,2]. It explains the absence of quarks 
from the physical spectrum by an analog to the MeiBner-Ochsenfeld effect of solid state physics: 
The vacuum contains condensed chromomagnetic monopoles^ which squeeze the chromoelectric 
field lines between two quarks into a flux tube. These vortex like structures can be identified best 
after reducing the gauge symmetry from the non- Abelian SU{N) to the Abehan U{l)^-\ Several 
such partial gauge fixings are known [5]. The most prominent one is the maximally Abelian gauge 
(MAG), which is defined by minimizing the components of the gauge field which are not Abelian. 
Since the corresponding generators of the algebra are off-diagonal matrices, these field components 
are also called off-diagonal gluon fields in contrast to the diagonal ones, which correspond to the 
Abelian part of the gluons. 

The dual superconductor picture suggests that for the infrared (IR) properties of Yang-Mills 
theory the diagonal gluons play a dominant role, since the monopoles live in the Cartan subalgebra 
[6]. This idea was formulated by Ezawa and Iwazaki as the hypothesis of Abelian IR dominance. 
In their original paper [6] the exact mechanism of this Abelian dominance was not specified any 
further. In terms of the Green functions one possible realization is a massive behavior for all 
three propagators (diagonal and off-diagonal gluons and ghosts), but with the lowest mass for the 
diagonal gluon such that the other fields effectively decouple and the IR dynamics is dominated 
by the diagonal field. Such a behavior was found by lattice calculations [7, 8] and in the refined 
Gribov-Zwanziger framework [9, 10]. Another possible realization of Abelian IR dominance is a 
diverging diagonal gluon propagator. Such a behavior was found in an IR scaling analysis of the 
MAG [11, 12]. We want to stress that the coexistence of both solutions is not unlikely. Indeed, 
such a scenario is also known from the Coulomb [13, 14] and the Landau gauges [15, 16]. 

We will recall the definition of the MAG in Section 2 and recapitulate the cornerstones of an 
IR scaling analysis in Section 3, where also the results for the MAG are presented. Additionally we 
provide new numerical data corroborating the possible existence of the scaling solution in the MAG 
and explain how the found scaling solution might be related to the results from lattice calculations. 
In Section 4 we make a short detour to linear covariant and ghost anti-ghost symmetric gauges. 

2. The maximally Abelian gauge 

For the MAG we spht the gauge field into diagonal and off-diagonal components as follows: 

A^ = T'A'^ + T"Bl, (2.1) 

where the T' are the A'^ — 1 generators of the Cartan subalgebra of SU (N), i.e., [T' ,7^] = 0. For 
example, in SU (3) the diagonal generators are and T^. A widespread convention is to use the 

'Although magnetic monopoles are gauge independent objects [3] their detection in lattice calculations depends on 
the choice of gauge [4]. 
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indices /, j and so on for diagonal values only and a, b and so on for off-diagonal ones. Furthermore 
we use A for the diagonal and B for the off-diagonal gluon fields. This splitting has interesting 
consequences for the possible interactions, which can be inferred from the standard commutation 
relation of the generators: 

[T\r] = if"r, (2.2) 

where r, s and t can be either diagonal or off-diagonal. From eq. (2.2) we can directly see that only 
three off-diagonal fields or two off-diagonal and one diagonal fields can interact. Even more, as in 
SU (2) there are only two off-diagonal generators the first possibility also drops out. Thus we have 
a higher number of interactions in SU{N > 2) than in SU{2). In summary the pure Yang-Mills 
part has the following interactions: ABB, AABB, BBBB, (BBB, ABBB). The two interactions in 
parentheses vanish for SU (2). 

The idea of the MAG is to minimize the norm of the off-diagonal gluon fields so as to stress 
the role of the diagonal part. The condition 

RMAG = \j dxBl{x)Bl{x) ^ min (2.3) 

is fulfilled if D°^B^^ = 0, where the respective covariant derivative contains only the diagonal gluon 
field: 

Df:=d"'d^+gf"<A'^. (2.4) 

Since this gauge fixing condition is non-linear the resulting Faddeev-Popov terms in the action 
contain quartic terms. The remaining U{1)^^^ symmetry of the action is fixed to the Landau 
gauge, d^A'^ = 0. Since the field A' is Abelian the corresponding ghosts decouple so that only the 
ghosts of the non-diagonal sector remain. Finally, a quartic ghost term is added to the action for 
renormalizability [17]. The interactions stemming from the gauge fixing and the renormalizability 
requirement are Acc, AAcc, BBcc, cccc, (Bcc, ABcc). The two interactions in parentheses vanish 
for SU{2). 

The complete resulting action is 

Smag = J ^^(^f^vf^v + \F;:,F;:,+c"DfD';[c'-gf^''c"DfB'=^c'^^ 

+ lg'af''f'-c"c'c'c' + a\g'f''f"'^c'(^c'c' + ^{d^A'^f) , (2-5) 

where a is the gauge fixing parameter of the non-diagonal part and ^ that of the diagonal one. 
Since a appears in the denominator of vertices it cannot, in contrast to the Landau gauge, be set to 
zero right away. The dependence on this parameter will be further elaborated in Section 3. 

Since the MAG has thi^ee fields with eleven interactions the derivation of its Dyson-Schwinger 
equations (DSEs) or functional renormalization group equations (FRGEs) is very cumbersome 
when doing manually. To this end the program DoFun [18] has been developed which can ex- 
actly do this for us and the derivation process is considerably shortened. In fig. 1 the DSE of the 
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Figure 1: DSE of the diagonal gluon two-point function obtained with DoFun [18, 19]. Continuous lines 
denote diagonal gluons, dashed ones off-diagonal gluons and dotted ones ghosts. 



gluon two-point function is shown as an example. DoFun is the successor of DoDSE [19]. Two 
new features are the inclusion of FRGEs and the derivation of the Feynman rules from a given 
action. Furthermore the complete algebraic form of the integrands can now be obtained so that 
computations directly with Mathematica are possible. This was very helpful for the calculations 
of the IR leading diagrams of the IVIAG [12] and the Gribov-Zwanziger action [20]. For the latter 
the application of DoFun helped to identify a unique solution [21]: Only with the help of a com- 
puter algebra system we could show that one of the two possible scaling solutions does not yield a 
numerical solution for the infrared exponent and can thus not be realized. 

3. Scaling analysis 

Since functional equations are formulated in terms of non-perturbative Green functions they 
are well suited for the investigation of the IR regime of Yang-Mills theories [22, 23]. The gauge 
which received the most attention is without doubt the Landau gauge. Here two solutions have 
emerged from the analysis of the system of gluon and ghost propagator DSEs: The scaling solution 
[16, 20, 22-28] and the decoupling solution [15, 16, 27, 29, 30]. The former is characterized by 
an IR vanishing gluon propagator and an IR enhanced ghost propagator, whereas the latter has an 
IR finite gluon propagator and a finite ghost dressing function. Both solutions are connected by the 
choice of a boundary condition required to solve the system of equations [15, 16]. Whereas such 
a choice might seem arbitrary at first sight, it turns out that it only affects the deep IR behavior of 
the propagators, while all solutions agree for higher momenta. Up to now no calculation exists that 
shows a difference between these two solutions for physical quantities: In ref. [31] it was shown 
that both solutions yield a confining Polyakov loop potential, the confinement-deconfinement and 
the chiral transition temperatures were found to be the same for both solutions [32] and even meson 
masses are unaffected by choosing one of the two [33]. Consequently, the interpretation of this 
boundary condition as an additional gauge fixing parameter as in the Landau-S gauges of ref. [34] 
seems reasonable. We will see that for the MAG we also have two qualitatively different solutions 
and a similar connection between them exists. 
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For the scaling analysis we make the ansatz that all dressing functions are characterized by a 
power law in the IR, e.g., for a propagator we have D{p) = c{p^) j and c'^ = d ■ (p^)^, where 5 
is the so-called infrared exponent (IRE). As we are only interested in the qualitative behavior of the 
Green functions we focus on the different IREs. Indeed, it is even possible to shift the analysis to the 
level of the IREs, because all integrals are dominated by low momenta if the external momenta are 
low and we can replace all propagators and vertices by the corresponding IR expressions [28, 35]. 
The scaling of such an expression can then be determined by counting all exponents of momenta. 
Thus we can determine the IRE of any given diagram. 

The best method to derive the scaling relations of the IREs is the combined use of DSEs and 
FRGEs as first done for the Landau gauge in refs. [36, 37]. During the study of the MAG it turned 
out that quite generically the infinitely many relations can be reduced to a rather small number. 
For details we refer to refs. [11, 12] where these general relations are written down. Knowing 
them even allows to investigate complex actions Uke the Gribov-Zwanziger action with functional 
methods [20]. 

For the MAG we parametrize the propagators as (i^ =0) 

D>^) = 5'V^(,,.-^), (3.1) 
Dfip') = 5"'^'-^ - (1 - a)^) , (3.2) 

P 

and take the following power laws for the dressing functions in the IR: 

ca(/) ' = 'd, . (/)^^ csip') '^^'ds ■ {p')'\ cAp') ' = "d,-{p'f. (3.4) 

As we do not set the gauge fixing parameter of the off-diagonal part to zero the longitudinal part of 
the off-diagonal propagator could have its own dressing function. In order to take that into account 
one could introduce an extra IRE dsjong and split the off-diagonal field in transverse and longitudi- 
nal parts. However, inserting the bare vertices explicitly and projecting the DSE of the off-diagonal 
propagator transversely and longitudinally one finds that the equations for both parts are the same. 
Thus the IR analysis yields 8b = SB.iong- This result can only be invaUdated if cancellations between 
different diagrams occurred or due to contributions from the unknown dressed vertices. Based on 
this relation we adopt in the following only one common dressing function for both tensors. 
The only consistent scaling relation which can be obtained for the MAG is [11] 

i<mag-=-Sa = Sb = 5c>0. (3.5) 

This means that the diagonal degrees of freedom are IR enhanced and the off-diagonal ones ai^e 
IR suppressed. An upper bound for the parameter Kmag can be obtained from demanding that the 
Fourier transformations of the propagators should exist [26] : Kmag < 1 ■ The IREs for the vertices 
can also be obtained [11]. It turns out that vertices become more IR divergent when the number of 
off-diagonal legs increases. The method developed in [1 1] has been applied to the MAG in SU (2) 
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and SU {N > 2). Interestingly, no differences have been found. The additional interaction of the 
latter do neither invalidate the SU (2) solution nor do they allow an additional one. 

In the derivation of this scaling solution the bare two-point functions play a special role in the 
respective DSEs. In order to allow for a scaling solution at least one of them has to drop out at zero 
momentum. This is well known for the Landau gauge, where the bare ghost two-point function is 
cancelled when choosing the corresponding renormalization condition [16, 24, 38]. As we do not 
know at the beginning which function could be IR divergent, we allow all of them to become IR 
divergent. At the end it turns out that the only consistent solution is the one with an IR enhanced 
diagonal gluon propagator. This requirement immediately opens the possibility to check whether 
the value of the diagonal two-point function at zero momentum can serve as an additional gauge 
fixing parameter similar to the ghost two-point function in the Landau -B gauges [34]. Indeed, 
when we take an IR finite diagonal gluon propagator we find that the off-diagonal propagators 
become IR finite too. The reason is the presence of quartic interactions between diagonal and off- 
diagonal fields of the MAG. If the diagonal gluon field is massive, the tadpole diagrams in the DSEs 
originating from these interactions are proportional to the diagonal gluon mass. Consequently 
the left-hand side should be constant for vanishing momentum from which we conclude that the 
off-diagonal propagators also become IR finite. This is the mechanism how the zero momentum 
value of the diagonal gluon propagator switches between the scaling and the decoupling types of 
solutions. Note that there is no quartic interaction between the ghost and the gluon in the Landau 
gauge and consequently an IR finite gluon propagator does not entail an IR finite ghost propagator. 

The result of the scaling analysis described above is of course not sufficient evidence that 
this solution really exists. A full-fledged numerical calculation, on the other hand, is complicated 
by the structure of the scaling solution of the MAG: The IR leading diagrams are all two-loop. 
Consequently, for a consistent numerical treatment a truncation scheme rather different from the 
one used in the Landau gauge has to be found. In the latter it is possible to neglect the two-loop 
diagrams and still have a truncation consistent in the IR and the UV, because the leading diagrams 
in both regimes are one-loop. 

As a preliminary step towards a numerical (scaling) solution in the MAG we calculated a 
value for the parameter Kmag- For this one only needs the IR leading diagrams. Here we find again 
an ambiguity, since it is not clear whether the squint diagrams are IR leading [11, 12]. We only 
took the sunset diagrams, which are definitely IR leading, so that the properly projected equations 
reduce in the IR to 

' = -^aabb{p^^ KMAG)dAdl-X^^Xp^, KMAG)dAdl, (3.6) 
ds ' = —XaABb{P^^ ^MAG)dAdB, (3.7) 

d~^ = —X^^^{p^,KMAG)dAdc- (3.8) 

The X{p^, Kmag) denote the sunset integrals without the coefficients from the propagator power 
laws. The superscript indicates the corresponding DSE and the subscript the bare vertex contained 
in the diagram. Using the invariant combinations Ii := d^dg and I2 '■= d\dl the three equations can 
be combined to 

XLbb(p^^^mag) X^^.^{p^,Kmag)' 
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Figure 2: On the left-hand side the evaluation of eq. (3.9) is shown for several values of the gauge fixing 
parameter a. Crossings with the horizontal line represent solutions for Kmag- On the right-hand side the two 
solution branches for varying a are shown. 

If the X{p^, Kmag) are known, this equation yields the solution(s) for Kmag- The computation of 
the integrals requires the dressed four-point functions. Their power law behavior is known (they 
are IR constant), but their tensor structures are not. As simple ansatze we employed the tree-level 
structures^. The last unspecified quantity is the gauge fixing parameter of the off-diagonal part, a. 
In the left plot of fig. 2 we show eq. (3.9) for several values of a. Surprisingly there are solutions 
0.7 < Kmag < 0-8 for all tested values of a. The plot on the right-hand side underlines this: For 
certain ranges of a there are two solution branches, but there is always one around Kmag ^ 0.74. 
The stability with respect to changes of a is by no means required and even persists for rather 
large values of a. This is in big contrast to linear covariant gauges, where the existence of a naive 
scaling solution in accordance with the Slavnov-Taylor identities is only assured for ^ = 0, see 
Section 4. We want to stress the fact that finding a sensible solution for Kmag is a non-trivial result: 
The scaling analysis itself does not rely on any details of the vertices except for their combinatorial 
properties, while the calculation of Kmag depends on all the little details such as Lorentz and color 
structure. In this calculation the program DoFun [19, 18] was of great help, since the intermediate 
expressions of the calculation would cover several pages. 

4. Other gauges 

The scaling analysis can be also be employed to other gauges. Here we present the results for 
linear covariant and ghost anti-ghost symmetric gauges, which generalize the analysis of ref. [39], 
where bare vertices were employed, to the case of dressed vertices. 

The decisive difference between the MAG and the linear covariant gauges is that the longitu- 
dinal part of the gluon propagator does not acquire a dressing in the latter but stays proportional to 
the gauge fixing parameter f : 

0?..(/)=«"(.„.-^f )^+«^^. (4.1) 

^All tree-level expressions can be found in [12]. 
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This can be inferred either from a Slavnov-Taylor identity [22] or from the fact that the longitudinal 
part of the gluon propagator, {{d^A^){dvAv)) , corresponds to the second moment of the Gaussian 
distribution, e 2f ("^f"**!) ^ ^j^^j ^gnce proportional to t, [40]. To include the longitudinal part of 
the gluon into the present formalism we can interpret it as an additional field. Its propagator has 
the IRE SA^iong- This allows to distinguish also between the transverse and longitudinal parts of 
the vertices. Naively the same inequalities are obtained for the longitudinal gluon field as for 
the transverse one. This could change if some vertices vanished when they are contracted with a 
longitudinal gluon propagator. After a more careful analysis one finds that indeed some diagrams 
drop out, but at the end those turn out not to be important. The vertex which is decisive for the 
IR leading behaviour is again the ghost-gluon vertex and we obtain for the longitudinal gluon 
propagator IRE the same result as for the transverse one, as the longitudinal part of the ghost-gluon 
vertex does not vanish. The result 8a = ^A.iong = —25^ leads together with the triviality of the 
longitudinal part, SA.iong = 0> to 5a = —25c = [11, 12]. This allows two possible conclusions: 
Either there is no non-trivial scaling relation for linear covariant gauges, or the naive application of 
the IR analysis is insufficient here. 

Finally we consider ghost anti-ghost symmetric gauges. Their most important difference to 
the Landau gauge is that they feature an additional quartic ghost interaction. This interaction is 
exactly what denies a consistent non-trivial scaling solution for this gauge [11, 12]: As in linear 
covariant gauges we obtain 8a = 8c = 0, i.e., a trivial scaling solution. 

5. Summary 

The MAG is a convenient choice of gauge for the investigation of the hypothesis of Abelian IR 
dominance. In addition to the decoupling type solution already known from lattice simulations and 
the refined Gribov-Zwanziger framework our results hint at the existence of a scaling type solution 
where the diagonal gluon propagator is IR dominant and the off-diagonal degrees of freedom are 
IR suppressed. The IREs for the complete tower of Green functions can be derived in this case. The 
two types of solutions are related via the zero momentum value of the diagonal gluon two-point 
function: If it is chosen as zero the scaling solution is obtained, whereas we could show that a 
finite value at zero momentum automatically leads to finite values at zero momentum for the off- 
diagonal propagators, i.e., the decoupling type solution is realized. If the situation is similar to the 
Landau gauge, this choice is a choice of gauge and should not have any consequences for physical 
observables. A common feature of both solutions is that the diagonal degrees of freedom dominate 
in the IR in support of the Abelian IR dominance hypothesis. Thus further investigations using 
the MAG are promising to obtain a better understanding of the confinement mechanism. Finally 
we want to stress that the MAG is the first gauge besides the Landau gauge where a full scaling 
solution was found. We investigated also linear- covariant and ghost anti-ghost symmetric gauges, 
but straightforward scaling solutions do not seem possible and more detailed analyses are required. 
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